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Abstract. Vignale and Kohn have recently formulated a local density approximation to 
the time-dependent linear response of an inhomogeneous electron system in terms of a 
vector potential for exchange and correlation. The vector potential depends on the induced 
current density through spectral kernels to be evaluated on the homogeneous electron 
gas. After a brief review of their theory, the case of inhomogeneous Bose superfluids is 
' considered, with main focus on dynamic Kohn-Sham equations for the condensate in the 

r— 5 linear response regime and on quantal generalized hydrodynamic equations in the weak 

inhomogeneity limit. We also present the results of calculations of the exchange- 
correlation spectra in both electron and superfluid boson systems. 
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1. Introduction 



o 

<3 The evaluation of the dynamical properties of inhomogeneous quantal systems is a cen- 

00 tral problem in condensed matter physics. As an example we may recall the current ef- 

forts at elucidating the dynamical behaviour of dilute alkali vapours which have under- 
gone Bose-Einstein condensation in magnetic traps. The experimental studies have con- 
cerned (i) the excitation of low-lying shape deformation modes in a regime where there is 
£h no detectable noncondensate fraction [1,2] and the behaviour of mode frequencies and re- 

I laxation times as temperature is raised towards and above condensation [3]; (ii) the prop- 

""O agation of sound waves in the condensate and its thermal cloud [4, 5]; and (iii) antisym- 

G metric oscillations of the condensate and the thermal cloud corresponding to second 

sound [5]. These experiments have stimulated a number of theoretical studies (see e.g. 
[6] and references given therein). The relationship between the dynamics of such confined 
fluids and that of a homogeneous Bose superfluid [7, 8] also is of great interest. 

Local density approximations have been very useful in accounting for exchange and 
correlation (xc) in the ground state energy of inhomogeneous electron systems, and the 
search for dynamic extensions has attracted interest for some time (see e.g. [9]). A scalar 
time-dependent xc potential, which is a local functional of the time-dependent particle 
density, is useful in dealing with low-frequency phenomena. However, analysis of the 
constraints coming from basic conservation laws has shown that inconsistencies can 
arise and are associated with the non-existence of a gradient expansion at finite frequency 
CO for the xc potential in terms of the density alone (see e.g. [10]). Recently Vignale and 
Kohn [11, 12] have overcome these difficulties by resorting to a dynamic xc vector po- 
tential even in the case of a system subject to an external scalar potential. They got an 



expression for the xc vector potential in the linear response regime in terms of kernels to 
be evaluated on the homogeneous electron gas at the local equilibrium density. Their 
expression becomes exact when the equilibrium density and the external potential are 
slowly varying in space, on length scales set by l/k F and v F /co where k F and v F are the 
local Fermi wave-number and velocity. The kernels are directly related to current-current 
response functions [13]. 

In § 2 below we give a brief review of the work of Vignale and Kohn and then pro- 
ceed in § 3 to report on our work regarding an inhomogeneous superfluid. Again for the 
latter the validity of the results is restricted to weak inhomogeneity (slow variations over 
length scales set by the interatomic distance and by c/co where c is the local speed of 
sound) in the linear response regime. In § 4 we report some examples of xc spectra for 
both the electron fluid and the Bose superfluid, which have been obtained in basically the 
same decoupling scheme for the current-current correlations [13]. 

2. Weakly Inhomogeneous Electron Fluid in the Normal State 

According to the Runge-Gross theorem [10] the problem of a many-electron system in 
an external potential V (r, t ) can be mapped into that of non-interacting electrons in an 
effective potential V(r,t) + V ex (r,t) (the excess potential V ex (r,t) includes both the 
Hartree term and the xc term). In the so-called "adiabatic local density approximation" 
(ALDA) the dynamic xc potential is evaluated from the xc energy density of the homo- 
geneous electron gas in the form 

V x WA {r,t)=de h xc {n{r,t))l dn , (1) 
as in the corresponding static problem. The ALDA embodies the correct electron-gas 
compressibility, but does not account correctly for plasmon dipersion and omits plasmon 
damping as well as transverse-current fluctuations. In the context of linearized hydrody- 
namics for a monatomic fluid, the ALDA embodies the speed of isothermal first sound, 
but does not account for sound-wave damping nor does it allow for the transition from a 
collisional to a collisionless regime which may take place e.g. as temperature is varied. 

As already remarked in § 1, Vignale and Kohn [11, 12] have proposed using a 
dynamic xc vector potential in the linear-response regime to overcome the limitations of 
the ALDA for the dynamics of an inhomogeneous electron system in the normal state. 
The induced current density is written as 

5j(r,ffl) = |rfi J x JfS (r,r';»)[a(r,ft)) + a, / (r,») + a xc (r,w)] , (2) 
where % KS (r, r' ; (o) is the Kohn-Sham response of a non-interacting reference system, 
a(r, co) is the applied potential, a H (r, co) is the Hartree mean-field term and a jx (r, co) is 
the xc vector potential. This is related to the induced current density by 

aJr,co) = \dr'fJr,r';co)dj(v',co) . (3) 

The kernels f xc (r, r" ;co) entering Eq. (3) have some general exact properties, which are 
needed to relate them in the weak-inhomogeneity limit to analogous kernels for the ho- 
mogeneous electron fluid taken at the local equilibrium electron density n(r) [12]. These 
are Onsager's reciprocity relation, the zero-force and zero-torque theorems (in brief, the 
net force and the net torque exerted by the xc potential on the system must vanish) and a 
Ward identity embodying the compressibility sum rule for the homogeneous electron 
fluid. As the final result the dynamic correction to the ALDA xc vector potential is 



expressed in terms of the induced velocity u(r, co) = <5j(r, co) I n(r) and of the homoge- 

neous xc current-density kernel f xc (k; ft)), evaluated at the local equilibrium density. The 
latter may be separated into the sum of a longitudinal and a transverse part: 

f^(k;ffl) = ffl" 2 [/^(ffl)kk + /^(ffl)(A; 2 I-kk)] . (4) 
Rather than reporting the complicated expression of the xc vector potential derived by 
Vignale and Kohn [11] in the local density approximation, we refer to the later work of 
Vignale et al. [14], who rewrite it in a physically transparent way through a visco-elastic 
stress tensor a . Defining 

<> = nJ co)(^ + ^L-^8^-n) + U^S,V-n (5) 

where t] xc ((0) and £ vt .(ft>) are complex viscosity coefficients given by 

riJco;n) = -(n 2 1 ico)f T Jco\n) (6) 

and 

C(6);n)-(» ! /ifi))[/>;»)-|/I(ffl;'>)-'' ! e,W/tf] - (7) 

one finds that the divergence of cP' 1 determines the dynamical xc vector potential. The 
connection of the theory with generalized hydrodynamics in the linear-response, long- 
wavelength regime becomes fully evident when one introduces a weak-inhomogeneity 
assumption also for the Kohn-Sham response matrix, by replacing it by its diamagnetic 
part Xks (r, r' ; <») -> <5(r - r' )n(r, co)/ m. The induced current is then given by 

^ z 
-imcodj(r,co) = n(r)[-icoa(r, co)- V(5v H (r, CO)- V(8n(r, ft)) / n 2 K r )] + V- c(r, co) (8) 

where 8v H (r, co) is the Hartree potential associated with the long-range Coulomb inter- 
actions, K T is the local electron-gas compressibility and the full stress tensor 6 differs 
from the expression in Eq. (5) by the replacement r\ xc (co) -> T\ x (co) -p (n)/ ico , p be- 
ing the pressure of the non-interacting electron gas. 

In summary, the corrections to the ALDA restoring force in the linear-response, 
weak-inhomogeneity regime merely involve the visco-elastic effects associated with fre- 
quency-dependent elastic constants and damping coefficients, which are related to the 
imaginary part and the real part of the two viscosity coefficients T] xc (co) and £ tt (ft)), re- 
spectively. For longitudinal motions, in particular, the plasmon excitation arises from 
the Hartree term in Eq. (8) and the role of the visco-elastic function £ vc (ft)) + 4rj xc (C0) 1 3 
is (i) to shift the plasmon dispersion coefficient away from the low-frequency value de- 
termined by K T in Eq. (8), and (ii) to introduce plasmon damping from decay into mul- 
tiple electron-hole pairs. 

Equation (8) relates the generalized viscosity coefficients ?7„(ft)) and £ xc (Cd), or 
equivalently the dynamic local field factors f xc (C0) and f T xc (co) entering Eqs (6) and (7), 
to the current-current response function of the homogeneous electron gas, via a finite-fre- 
quency generalization of the Kubo relations [13]. This property will be used in § 4 for an 
evaluation of these xc spectra, based on the exact equation of motion for the current re- 
sponse and on the reduction to two-pair excitation processes through a decoupling ap- 
proximation [13, 15]. We shall first turn, however, to a discussion of an inhomogeneous 
Bose superfluid. 



3. Dynamic Density Functional Theory for a Bose Superfluid 



In a superfluid of Bose particles the set of basic dynamic variables must include the 
superfluid velocity field u,(r,f) in addition to the total (superfluid plus normal-fluid) 
current density. We follow Hohenberg and Martin [7] in introducing the transformation 

Y i (r) = [n (r)]" 2 exp[-/<p(r)] (9) 
for the field operator in terms of the condensate density and phase operators: then the su- 

J_ 

perfluid velocity operator is u s (r) = m~'V<p(r), in accord with the irrotational character 
taken by the superfluid velocity field in the absence of vortices. 

The linear response of the condensate density and phase to a dynamic gauge -breaking 
external field can be treated by means of the equation of motion for the field operator and 
is found to have the Hohenberg-Kohn-Sham structure compatible with a mapping of the 
interacting system into a single-particle reference system [16]. The condensate self- 
energy a(l) is introduced by 

C7(l) = jrf2v(r 1 ,r-)<^ t (2)^(2)^(l)> (10) 
where l = (r 1 ,t 1 ), v(r 1 , r_ ) is the interparticle pair potential and the symbol <•••> de- 
notes the expectation value of an operator on the equilibrium ensemble at given tempera- 
ture. The Kohn-Sham response functions of the condensate are constructed 
from a reference Green's function G KS (1, 1' ) which is given by 

G^(U') = [/-f- + -^V 2 -V(r 1 ) + /i g ; (ri) ]S(l,Y) , (11) 

where V(r,) is the confining potential and fi is the chemical potential. The two-by-two 
matrix %{\, 1' ) of the response functions of the condensate can then be cast into the form 

X = X ks +Xks® k ®X (12) 
where the symbol ® denotes integration over intermediate variables and the kernels in 
the matrix K are obtained as functional derivatives of the real and imaginary parts of 

Sa(l)=a(l)-a eq (v 1 )< ¥ (l)>/< ¥ (r l )>^ ~ (13) 

with respect to the condensate density and phase. 

We should emphasize that (i) the Kohn-Sham reference system that we are proposing 
for the superfluid contains the interactions through the equilibrium value of the conden- 
sate self-energy entering the RHS of Eq. (11); and (ii) we are not advocating a density 
functional approach to the evaluation of the equilibrium state of the superfluid. In partic- 
ular, a thermodynamic treatment based on the Gross-Pitaevskii equation already is in 
quantitative agreement with experiment for confined vapours of alkali atoms which have 
undergone Bose-Einstein condensation [17]. 

A similar fully microscopic treatment of the current density induced in the inhomo- 
geneous superfluid by a further external driving field leads to microscopic expressions for 
the remaining response functions (the current-current response and the current-condensate 
cross response) [16]. While these expressions are quite complicated in the general case, 
these calculations allow one to complete contact with the well known two-fluid model of 
Tisza and Landau. Firstly, the total current-current response function has the general 
structure already found by Hohenberg and Martin [7], who used it to show that in an 
appropriate limit this function gives the response of the total current density to an 
external field defining the normal fluid velocity u n (r, t) . Secondly, a Ward identity for 



the current response of the noncondensate at fixed condensate, which was first given by 
Huang and Klein [18] can be used to introduce a function p,(l,2) through the equation 

V 1 p,(l,2)=V l [n t .(l),5(l,2)]-2mIm[<^ t (l)>|^l < , /> ] . (14) 

oa(2) 

Within our microscopic treatment Eq. (14) can be rewritten in the form 

Vj> (1,2) = V . (15) 

11 2 <5u s (2) " 

It follows from Eq. (15) that the function p,(l,2) reduces for the homogeneous fluid in 

the static limit to the hydrodynamic definition of superfluid density given by Hohenberg 

and Martin [7]. It can also be shown that in the same limit the definition of superfluid 

density given by Griffin [8] is recovered from Eq. (14). 

On the above grounds we have developed a two-fluid model for the generalized hydro- 
dynamics of the inhomogeneous superfluid in the linear-response, weak-inhomogeneity 
regime [16]. We propose that Eq. (14) may provide a reasonable definition of the equilib- 
rium superfluid density p s (r) in the weakly inhomogeneous case, when we take its 
static limit and the k = component of its Fourier transform with respect to r 1 —r 2 . 
The functional derivative entering the RHS of Eq. (14) is a five-point correlation func- 
tion, to be evaluated in this case on the equilibrium state. 

With the above definition of the superfluid density, we first proceed to extend 
Landau's hydrodynamic theory for the homogeneous superfluid to finite-frequency phe- 
nomena. We use for this purpose the well-known memory function formalism, as de- 
scribed e.g. in the book of Forster [19]. We assume isothermal conditions, i.e. we 
neglect the couplings between temperature and density fluctuations. The form of the gen- 
eralized hydrodynamic equations is dictated by some general considerations: (i) invariance 
under a Galileian transformation and Onsager symmetry must hold; (ii) as a consequence 
of the zero-force and zero-torque theorems, the time derivative of the current density j is 
driven by the divergence of a symmetric tensor of the second rank; (iii) the time deriva- 
tive of the superfluid velocity u is the gradient of a scalar, in view of its irrotational 
character below threshold for vortex generation; and (iv) the internal driving forces are 
determined by the divergence of the normal-fluid velocity u n and of the interdiffusion 
current j r = j - nu n = p, (u t - u n ) . 

The generalized hydrodynamic approach is easily extended to a weakly inhomoge- 
neous superfluid in isothermal conditions [16]. The general properties that we have listed 
above to infer the structure of the generalized hydrodynamic equations in the homoge- 
neous case remain valid. In addition, a Ward identity is essential in relating the effects of 
a weak inhomogeneity on the excess kernels to their density dependence: it allows the 
inhomogeneous kernels to be set equal to those of the homogeneous fluid at the local 

densities of superfluid and normal fluid. The resulting dynamical equations are 

^ z 
-imcodj(r,co) = n(r)[-z'£oa(r, co) - V '(8n(r, co) I n 2 K T )] + V- c(r, co) (16) 

and 

-imcoSn s (r, co) = -icoa s (r, co) 

-Wn 2 K T )- > dn(r,co)-(Ts/c v )ds(r,co)] + V-(7< s) (r,CO) , (17) 
where s and c v are the entropy and heat capacity per particle, and the stress tensors are 
given by 



° v =[ri{(o)-p {n)l icot^ + " f <5,V • u„ ) 

dr or 3 

+ 5 iJ [f 2 (ffl)V-u„ + f 1 (ffl)V-jJ (18) 

and 

ff« = <^ 3 (ffl)V.j r + C 4 (ffl)V.uJ , (19) 
/?„(«) being the ideal-gas pressure. The ALDA contributions have been explicitly writ- 
ten in Eqs. (16) and (17) in terms of the total particle density fluctuations and of the en- 
tropy fluctuations, the latter being associated with fluctuations in the density of the non- 
condensate. The corresponding thermodynamic coefficients are to be evaluated on the 
homogeneous fluid at the local equilibrium density. Evidently, the ALDA driving forces 
are responsible for the first and second sound modes. 

The visco-elastic coefficients in Eqs. (18) and (19) are related to a set of excess 
kernels according to the following equations: 

= -(i£»r [/•„„, {(0)-d Pxc l dn\ T ] , (20) 
CA(0) = -(icor , [f, LJL ((0)-4f iTjT (co)/3-ndpJdn\ T ] , (21) 
£(£») = -{i(oy\f M J(D)-dii^ I dn\ T -(Ts/c v )d!U, c I dT\ n ] , (22) 
{ 4 (co) = -(icoy[f ujL (co)-dpJdti\ T ] (23) 

and 

ri((o) = -{i(orf iTiT {(o) , (24) 
with = £ 4 from Onsager symmetry. As already remarked, the kernels in these equa- 
tions are those of the homogeneous fluid at the local equilibrium densities of superfluid 
and normal fluid at temperature T . The xc pressure and chemical potential are similarly 
those of the homogeneous fluid at the local equilibrium density. 

Within the memory function formalism, the visco-elastic spectra in Eqs. (18) and 
(19) can be calculated from the linear response functions of the homogeneous superfluid 
through finite-frequency forms of the Kubo relations. These will be used in § 4.2 for a 
Bose vapour. 



4. Evaluation of dynamic kernels 

4.1. THE ELECTRON FLUID IN THE NORMAL STATE 

As we have already remarked at the end of § 2, Eq. (8) implies that the visco-elastic coef- 
ficients are related to the current-current response function of the homogeneous electron 
fluid by a finite-frequency extension of the Kubo relations. Precisely, 

l m f^\co) = lim^^^lmx UT) (k,co) , (25) 
n k 

where % UT) (k, ft)) is the proper longitudinal (or transverse) current-current response [13]. 
This emphasizes that the plasmon does not contribute to the xc kernels. 

Starting from the equation of motion for the current density fluctuation operator jr. , 
the following exact long-wavelength result can be derived: 



Im X„ (k, CO) = -(m 2 (0 4 )" £ £ r , (q, k)r. . (q' , -k) 

q,q' «' 

x Im{?J" o " <fo exp[i(fi) + i7/)f] < j^Pj ] >} + o(k 2 ) , (26) 

p r being the density fluctuation operator and the detailed expression of the coefficients 

r being given in ref. [15]. Decoupling of the four-point correlation function entering 
the RHS of Eq. (26) into all possible products of two-point correlation functions leads to 
the following approximate result [13]: 

T run, , r d(0' t d 3 q (4ne 2 1 nco' ) 2 T 

Imf T > (C0) = - J — \ y — — !- \m% L (q, CO - CO' ) 

% {2%f (co -co) 

x[a L{T) lmx L (q,0)' ) + b Lm (0)' /0)) 2 lmx T (q,0)' )] . (27) 
Here, the a's and b's are numerical coefficients given by a L = 23 a T 1 16 = 23 / 30 and 
b L = Ab T 1 3 = 8 / 15. Equation (27) only includes the contribution from excitations of 
two correlated electron-hole pairs. It can also be derived by means of second-order pertur- 
bation theory followed by approximate resummation to infinite order through renormal- 
ization of single-particle excitations from bare to screened [20]. For further details of this 
calculation, including the various possible choices for the response functions in the RHS 
of Eq. (26) and the inclusion of final-state exchange processes restoring the exactly 

known coefficients of the co 3 ' 2 decay of the spectra at high frequency [21], the reader is 
referred to the original papers [13, 15]. 

Equation (27) shows that both the longitudinal and the transverse spectrum involve 
an LL and an LT channel, the former becoming dominant at high frequency and the latter 
at low frequency (on a frequency scale set by the plasma frequency CO pl ). The values of 

the numerical coefficients in Eq. (27) are such that Im/J t = (16 /23)Im/^ at high fre- 
quency and Im f T Kc =(3/4)Im f L xc at low frequency: as a consequence, the two spectra 
have very similar shapes, the relation Im f T xc = 0. 72 Im f L xc being approximately valid at 
all frequencies. This is shown in Figure 1, which reports the two spectra for the electron 
fluid at a value of the coupling strength r = 3 in comparison with an earlier attempt by 
Gross and Kohn [22] at estimating the longitudinal spectrum by a smooth interpolation 
between high and low frequency. In contrast, our microscopic calculations show a sharp 
spectral threshold at 2 CO which is due to the opening of a two-plasmon channel in the 

two-pair excitation spectrum [23]. The large spectral strength of the plasmon excitation 
at long wavelengths, as compared to single-pair excitations, accumulates most of the os- 
cillator strength for two-pair processes in the spectral region just above 2 (O pl . 

The real part of the longitudinal and transverse kernels is obtained from the spectra il- 
lustrated in Figure 1 by means of Kramers-Kronig relations, i.e. via Hilbert transform 
and fulfilment of the third-moment sum rules. As a consequence of the spectral threshold 
at 2cO pl , a sharp minimum is found in Re/^ m at this frequency. Strong xc effects are 

present in Re f^ T) in the region of the plasma frequency, where these functions decrease 
with increasing CO towards their minimum at 2 CO . . This theoretical result is in accord 

with the available data on the plasmon dispersion coefficient from Electron Energy Loss 
experiments on the alkali metals [24], although band-structure effects will also play a 
significant role in a fully quantitative comparison between theory and experiment. 
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Figure 1. \mf^f((o) for the electron gas at r = 3 in units of 2(0 pl I n, as functions of 
(O I C0 pl . The dashed line gives the Gross-Kohn interpolation scheme. 

As a last point we briefly record the main results that have been obtained by this ap- 
proach [15] for the elastic moduli and the viscosity coefficients of the fully degenerate 
electron gas in dimensionalities D = 3 and D = 2. These are: (i) the calculated electron- 
gas compressibility is in good agreement with the available Quantum Monte Carlo data 
on the xc energy; (ii) xc effects decrease the shear elastic constant for r > 6 in D = 3 and 
at all r in D = 2; (iii) the bulk viscosity vanishes identically at zero temperature; and 
(iv) the shear viscosity takes a finite value on inclusion of dynamic xc effects and this 
value rapidly decreases with increasing r . 

4.2. THE DILUTE BOSON SUPERFLUID 

The visco-elastic spectra generalizing sound wave attenuation and shear viscosity to 
finite frequency are related by Eqs. (16) and (18) to the longitudinal and transverse 
current-current response functions of the homogeneous fluid. Precisely, these relations 
are 

4 — (am 2 

Re[C 2 (o» + - ri(co)] = lim^„ —^-Jm Xjlll (*, ffl) (28) 

and 
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Figure 2. Spectra of shear viscosity (left scale) and bulk viscosity (right scale) as 

functions of (0 1 (O at T = (full curve) and T - l.5(0 o (dashes). The dots 
show the result obtained at T = without a cutoff in the potential. 



Re J 7(©) = lim^ -^Im^ rJr (fe,«) . (29) 

Of course, the imaginary parts of £ 2 (ft)) and J7(ft>) , which are finite-frequency bulk and 
shear elastic moduli, are related to the real parts by Kramers-Kronig relations. Equations 
(28) and (29) are evidently analogous to Eq. (25) and can be evaluated in the same 
approach which led to Eq. (27) for the electron fluid. For the boson superfluid this 
approach may also be viewed as amounting to an expansion of the single-particle self- 
energy to second order in perturbation theory, followed by renormalization of the single- 
particle propagators into Bogolubov phonons. The extension of the theory to finite 
temperature is easily achieved, by introducing in the calculation of the response 
functions the thermal occupation factors Z?(ft>) = [exp(/?ft>)- 1]~' as required by the 
fluctuation-dissipation theorem. 

The case of present interest concerns the confined condensates of Na or Rb atoms in- 
side magnetic traps, in which the atomic interactions can be modelled by a repulsive 
Fermi model potential involving an s-wave scattering length a as a single parameter. To 
regularize the interatomic potential v k at large momenta k, we introduce a cut-off by 
writing 

v t =(4^a/m)exp(-fe 2 / 2(7 2 ) . (29) 
The corresponding Bogolubov spectrum is 

C0 k =[( W( ,/fc 2 /m) + (fc 2 /2m) 2 ]" 2 . (30) 



A lengthy but straightforward calculation, whose details will be reported elsewhere [25], 
yields for the shear and bulk viscosity spectra in the collisionless regime the result 

_. r , 12 „ „ v\k 6 2 3 2 s sinh(/?ft),) 

Re[ T?(fl»] = — Re[C 2 (a»] = — ' , (I ~ ~ a k + — a\ ) ^ * 7 . (31) 
5 240otm <u t w t 7 35 cosh(pco ( ) - 1 

Here, k is defined by 2(O k = CO, (O k ' = d(O l I dk and OL k — (k- Vv t ) / v k . At low fre- 
quency we find ?7((w) oc <y 3 if /J<y»l and ri{(o)<*(Q 2 T if /J<w«l. 

Figure 2 shows these visco-elastic spectra as functions of frequency at T = and at a 
higher temperature on a scale set by the frequency co = \nna I m . A further result is 
obtained by evaluating the longitudinal spectrum at the sound-wave frequency ck , 
whereby one obtains the width A (k) of the phonon as 

A(k) = 3k 5 1(640 Ttmn) (32) 

at T = and 

A(k) = 37t'k(k B T) 4 /(40mc 4 ) (33) 
for pco«l. These results for the phonon linewidth are well known from earlier work of 
Belyaev and Popov [26]. 

It can be shown that the results reported above for the visco-elastic current spectra at 
T = correspond to the so-called one-loop approximation introduced in the context of 
Bose gases by Wong and Gould [27]. We have calculated the other spectra of relevance to 
Eqs. (16) - (19) by the same method, starting from the finite-frequency Kubo relations 

Re[?,(£0)] = lim^^lm^.. (*,£») (34) 
k 

and 

Re[C(ft))] = Re[C 4 (ft))] = lim^ — lm X , L „(k,co) . (35) 

It turns out that within the one-loop approximation all these dissipation spectra can be 
expressed in terms of four exchange-correlation building blocks: a proper condensate ker- 
nel related to the shift in local chemical potential due to interactions in the fluid away 
from equilibrium, two noncondensate kernels (the irreducible proper parts of the longitu- 
dinal and transverse current-current response) and a cross condensate-noncondensate vertex 
function. Actual calculation of these spectral functions at T = shows that in a colli- 
sionless regime they take the same value aside from simple multiplying factors. The 
results are 

yRe[C 2 («)] = |Re[m»C 1 (®)] = ^Re[m 2 n 2 C 3 («)] = Re[J7(»)] ■ (36) 



5. Summary and concluding remarks 

We have reviewed the progress recently made in developing the theory of the dynamics of 
weakly inhomogeneous quantal fluids in the linear response regime from knowledge of 
the dynamics of the corresponding homogeneous fluid. The present focus in this area is 
on the degenerate electron gas and on the dilute Bose superfluid, for which the basic theo- 
retical concepts and some inputs for applications have been developed. The fruitfulness 
of this approach remains to be tested through applications to specific physical problems. 
Such applications are already under way in the case of electron fluids [28] . 
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